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Abstract: A small set of final prototypes of the ATLAS Inner Detector silicon tracker (Pixel and 
SCT) were used to take data during the 2004 Combined Test Beam. Data were collected from runs 
with beams of different flavour (electrons, pions, muons and photons) with a momentum range of 2 
to 180 GeV/c. Four independent methods were used to align the silicon modules. The corrections 
obtained were validated using the known momenta of the beam particles and were shown to yield 
consistent results among the different alignment approaches. From the residual distributions, it is 
concluded that the precision attained in the alignment of the silicon modules is of the order of 5 /im 
in their most precise coordinate. 

Keywords: Detector alignment and calibration methods. Solid state detectors. Particle tracking 
detectors. Large detector systems for particle and astroparticle physics. 
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1. Introduction 

This note reports the results of the alignment of the ATLAS Inner Detector [1] silicon tracker (Pixel 
and SCT) modules at the ATLAS Combined Test Beam data-taking (CTB) which took place at the 
CERN H8 beam-test facility in 2004. The purpose of the CTB was to study the combined perfor- 
mance of ATLAS. The setup represented a full barrel slice of the Inner Detector (ID), Calorimeter 
and Muon Spectrometer of the complete ATLAS detector and was instrumented with final proto- 
types. 

Once the Pixel and SCT modules had been installed in the CTB setup in addition to the already 
operational TRT, the Inner Detector was fully integrated into the common data acquisition system. 
Data were collected with this fully integrated ID, using beams with different characteristics. Pion, 
electron, muon and photon beams were used in a wide range of momenta from 2 to 180 GeV /c, 
and some data were taken without magnetic field (B). 

The CTB setup represented an ideal framework for testing the Inner Detector software. The 
offline reconstruction was tested on real data using the ATLAS software framework (ATHENA) [2] 
and was particularly useful for tracking [3], pre-commissioning tests, and for testing the alignment 
software. 

Determining the locations of the tracking detector elements is crucial for the performance of 
the ID tracker. For this purpose, various alignment algorithms, based on optimization of track 
hit residuals, were applied to align the CTB silicon setup. An alignment algorithm specifically 
developed for the CTB (hereafter referred to as Valencia approach [4]) had been adapted from an 
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algorithm used in previous SCT standalone test beams [5], by the time the first data were collected. 
The Valencia approach produced alignment corrections for the initial CTB data analysis. For the 
final analysis of the alignment, three more algorithms were tested. These algorithms, developed 
for the alignment of the entire Inner Detector sihcon tracker, are: Robust [6], Local X'^ [7, 8] and 
Global [9, 10] approaches [11]. 

The resulting sets of alignment constants were used to measure the momenta of the incident 
particles in electron and pion runs. A comparison with the nominal momenta was used to cross- 
check the different alignment procedures. The residual distributions and reconstructed track pa- 
rameters were studied for electrons and pions with and without B field. The global reference frame 
was also studied by matching the alignment results via a global offset optimization. 

2. Setup, Data Samples and Tracking 

The Inner Detector volume in the CTB setup was divided into three containers for each sub- 
detector: Pixel, SCT and TRT. Six Pixel and eight SCT modules were placed in their respective 
containers' . The TRT setup consisted of two barrel wedges, equivalent to 1/16 of the circumference 
of a cylinder. 

The coordinate system was chosen to be right-handed, with the X-axis along the beam direc- 
tion and the F-axis pointing vertically upward as depicted in Fig. [I] [12]. The origin was located at 
the entrance of the dipole magnet that produced a maximum 1.4 T field in the negative Z-direction. 
The Pixel and SCT detectors were located inside the magnet whereas the TRT detector was located 
outside due to its larger dimension. 

Inner Detector Geometry 



-1000mm 




■1710mm- 



Figure 1. Schematic representation of the ID components and the magnet at the CTB. The reference coor- 
dinate system is also shown. The long pixel coordinate and the SCT module strips are almost parallel to the 
Z-axis. 

A Pixel module [13, 14] consists of a single silicon wafer with an array of 50x400 jxra^ pixels 
that are read out by 16 chips [15]. The active area of each module is ~ 60.8 x 16.4 mm^. In the 
CTB setup, six Pixel modules were distributed in three layers (0,1,2) and two sectors (0,1). The 
distances along the beam axis between the different layers and the locations of modules within 
each layer mimic the arrangement of the modules in ATLAS. The first Pixel layer was nominally 

'The ATLAS detector has, in total, 1744 Pixel modules and 4088 SCT modules. 
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located at 195.986 mm from the global coordinate center along the beamline (X-axis) and the last 
layer was located at 268.277 mm. Each module was positioned at an angle of about 20° with respect 
to the incident beam, around the long pixel coordinate. Modules in the same layer overlapped by 
~ 200 /im. 

A SCT module is built from four single-sided silicon microstrip sensors glued back to back in 
pairs with 40 mrad stereo angle for a 3D space-point reconstruction [16, 17]. The modules produce 
two hits, one in each plane. The SCT end-cap modules have a wedge-shaped geometry which 
results in variable pitch sizes (Fig. ^. In the CTB setup, one of the four shape-wise distinct SCT 
end-cap module types was used (outer module). For the outer end-cap modules, the readout strip 
pitch is 70.9-81.1 /xm. Each plane has a length of about 120.0 mm and bases of about 72.0 mm 
and 57.0 mm. The readout is provided by a binary chip [18]. Eight SCT modules were used in 
each of the four layers (0,3) of the CTB setup; distributed in two sectors (0,1) with a 4 mm overlap. 
The arrangement of the modules was similar to the SCT barrel configuration in ATLAS^, however, 
the modules were not mounted at an angle with respect to the beam axis. The SCT modules were 
nominally positioned from 378.198 mm to 598.218 mm along the beam axis. 

The beam-line instrumentation, including trigger and veto scintillators, Cherenkov counters 
and readout system is documented elsewhere [19,20]. The Inner Detector magnetic field profile 
was measured [12] and its non-uniformity was taken into account during the track reconstruction. 
The absolute momentum as measured by the silicon detector which was located in a very uniform 
magnetic field region was certified to better than 1 % by comparing the momentum reconstructed 
from silicon alone with that obtained independently using the angular measurement in the TRT. 

The CTB ID data taking was divided into five different periods between September 2004 and 
November 2004 [12], where 22 million usable events were collected. In order to evaluate the 
material effects in the tracker, aluminum plates (10% Xq) were inserted and removed between the 
Pixel, SCT and TRT setups (Fig. [ij) in alternate runs. The TRT was repositioned in the transverse 
plane of the beam. Particle type and energy of the beam also alternated during the periods. 

The algorithms provided a valid silicon detector alignment for all the CTB data-taking peri- 
ods. However, this article reports on the last period (period 5) of stable data-taking when no extra 
material layers were used. Table |l| lists the runs used for alignment studies in this period. Events 
from run 2102355, a 100 GeV pion beam run without a B-field, were used as input to all algorithms 
for the production of alignment corrections. For the Local approach, two other pion runs were 
used in addition. Further event selection details are given in Section ^. 

2.1 Simulation 

The CTB setup was simulated with Geant 4 using the same geometry description as the event recon- 
struction. Detector positions and initial numbers were provided through an Oracle-based conditions 
database (look-up information) which allowed the five different periods to be distinguished from 
one another. 

CTB specific modifications were applied to the simulation for studying the Pixel and SCT 
alignment, i.e. the propagation through material upstream of the ID and the inclusion of mea- 

^The rectangular barrel modules which have uniform 80 /in? pitch were not used due to their unavailability during 
test beam data-taking. 
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Table 1. List of selected runs used to assess the alignment results. 
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sured beam profiles. The upstream material (mainly air and triggering/monitoring scintillators) 
corresponded to 13.2% radiation lengths and was taken into account to mimic the momentum dis- 
tribution in the data properly. Profiles, consisting of beam incidence positions and angles, were 
taken from the data and were applied during the upstream simulation to bring the simulated hit 
maps and residual distributions of the silicon modules into agreement with the data. 

The magnetic field map was calculated taking into account the magnet geometry, in one quad- 
rant of the transverse plane with respect to the beam axis. The remaining field map was modeled 
assuming a symmetric the field map around the main axis of the magnet. The field map calcu- 
lated along these lines compares well with the actual measurement of the dipole field which were 
performed before and after the CTB runs. 

2.2 Tracking and Reconstruction 

The default tracking algorithm in the CTB was the 'CTBTracking' algorithm [3]. CTBTracking 
consists of a pattern recognition part, developed specially for the CTB, and a track fitting algorithm 
that is in use in full ATLAS as well as in the CTB. The pattern recognition finds the tracks by 
looping through combinations of space points. The track fitting algorithm is based on a global 
minimization technique, often called the 'breakpoint' method in the literature [21]. Multiple 
scattering and energy loss enter into the algorithm as additional fit parameters at a given number 
of scattering planes. The track fit has a custom description of the detector material in the test beam 
setup, with one scattering plane for each layer of siUcon modules. This material description was 
precisely tuned to give the best possible track resolutions, down to very low energies ( 1 GeV or 
less). A number of options and features exist in the track fit that are particularly useful for the 
alignment algorithms, such as the possibility of setting the momentum to a fixed value in the fit, 
and the ability to retrieve the fitted scattering angles and their covariances. 

Fig. ^shows typical hit maps for a Pixel and SCT module. The illumination was rather uniform 
for the channels that lay within the scintillator trigger acceptance window in the central region 
(~ 3 X 3 cm^ wide). More details on the tracking performance of the pixel detectors can be found 
elsewhere [22]. Unmasked noisy channels can be distinguished in the SCT bitmap. Those that were 
masked during data-acquisition appear as zero-entry channels. The illumination was not uniform 
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Col number Strip number 

Figure 2. Left: Pixel hit map for 100 GeV pion runs. The jc-axis corresponds to Pixel t] -coordinate and 
the y-axis to Pixel ^-coordinate. Empty horizontal bands correspond to the ganged pixels. Lighter vertical 
bands are due to the 600-^m-wide pixels. Right: SCT hit map for the same run. 



and limited along the strip length but only in the central region, where the sensor planes overlapped 
completely with the trigger scintillator. 

The limited illumination of the sensors had direct consequences on some of the alignment 
degrees of freedom (DoF) due to insufficient constraints and reduced sensitivity. The problem 
was more severe for SCT modules, because the SCT modules were not tilted with respect to the 
beamline. As the beam incidence was almost perpendicular to the module planes, the alignment 
procedures were not very sensitive towards misalignments along the beam axis. 

The pixel sensors require free space in order to bond the readout chips on the surface of the 
sensor. In the precise <p coordinate, unbonded pixels are physically connected to nearby pixels 
(ganged pixels) and share a readout logic channel. Due to this connection, whenever a hit was 
registered by a logic channel, there was an ambiguity as to which pixel fired. In the long coordinate 
wider pixels (600-/im instead of 400-/xm) are used. The wider pixels collect more hits. The impact 
of both effects is clearly seen in the pixel module hit map (Fig. |2|). The ambiguity in the ganged 
pixels was also found to effect the alignment. In a highly misaligned environment, tracking may 
make too many wrong decisions between ganged pixels. It was found that, in the presence of a high 
track quality cut, the ganged pixel hits were favoured, degrading the quality of the alignment [6]. 

The fact that the modules were exposed to almost perpendicular beams resulted in discrete 
Pixel Tj -residual distributions. Due to the large dimension in this direction (400 /im compared to 
300 jJ-m of the thickness of the silicon bulk) the drift of the charge carriers along that direction is 
negligible. Therefore, almost all of the clusters consist of a single pixel in the t] -coordinate. As 
the cluster position is located in its geometrical center, the outcome is a discrete positioning of 
clusters (Fig. ^). With only three pixel layers providing three precision points, a discrete residual 
distribution was obtained. The use of SCT clusters in the tracking partially removed this unde- 
sired effect [4]. Effectively the pixel tj -residuals of the first and last pixel layers were somewhat 
broadened by overlaps of Gaussian distributions, while the middle layer T] -residuals remained dis- 
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Figure 3. Pixel T] -residuals for tracks reconstructed with Pixels and SCT. Left: first Pixel layer. Right: 
middle Pixel layer 



Crete. This peculiarity of tlie CTB setup made the alignment along the pixel r\ -coordinate difficult. 
ATLAS collision data will not present such difficulties. 

3. Alignment of the CTB Data 

The goal of alignment is to determine the corrections to the parameters that describe the position 
and orientation of the module in space. Each module is treated as a flat rigid body with 6 DoFs, 
i.e., three translations along the local coordinate axes (x, y, z) and three rotations (a, j3, y) around 
the local coordinate axes, in a right-handed orthogonal frame where the origin is at the center- 
of-gravity of each module and the local x-coordinate is along the most precise coordinate. The 
translations correspond to the shift of the module with respect to its nominal position. For the axes 
orientation, the Cardano representation of angular rotation with respect to the cartesian axes was 
used. The alignment corrections were stored in the conditions database. 

The alignment corrections are given in terms of CLHEP [23] transform objects H, made of a 
rotation matrix R and a translation vector T . The rotation matrix is defined as: 

R = R,{a)-R,m-R-Sy) (3.1) 

with a, P and / being the rotation angles around the x, y and z-axes. Rz{Y) is the first rotation 
applied and Rx{oc) the last. The representation of a point in the local reference frame (P) of a 
module is H-P = f + RP in the global frame. Lets consider Hq as the transformation specifying 
the nominal position of a given module. If dH is a shift of the module, the new transformation of 
the points measured by it becomes H = Hq- 5H. Therefore, the task of the alignment is to determine 
the 6 DoFs that define 5H for each module. In the case of poorly constrained movements, some 
DoFs may not be considered. 

The technique to align each silicon module consists of minimizing its two residuals (pixel 
modules measure two coordinates and each SCT module has two sensor planes). The i-th residual 
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(defined by coordinate, plane or module) is thus r, = (m,- — e,(a,7r)) - ki, where represents the 
position of the hit recorded in the sensor plane, ^, is the intersection point of the extrapolated track 
to the detector that depends on the alignment parameters (a) and the vector of track parameters 
n = {do, 0o,Zo7 ^0) '?//')■ ki denotes the unit vector of the measurement direction [9]. 

All alignment algorithms were run iteratively. Initially, nominal detector and hit positions 
were used for track reconstruction. After the track fit, residuals and their derivatives with respect 
to alignment and/or track parameters were calculated to determine the alignment corrections. For 
each module, the best fit estimates for alignment parameters were derived and its position was 
updated. A new reconstruction with updated module positions was performed and the alignment 
was reiterated. This procedure is expected to converge to final alignment corrections for each 
module and the residual resolution is expected to improve. 

The alignment was performed using two different classes of approach. The Robust approach 
is based on iterative minimization of the residual means of overlapping and non-overlapping mod- 
ules. The approach is "robust" because the output is stable against changes in the input tracking 
information. 

The Valencia , Local and Global approaches are based on the linear least squares mini- 
mization defined for a set of reconstructed tracks as: 

xHa,7Cu...,nt)= £ Ti^Vr'n (3.2) 

/ € tracks 

where r,- = r, (a,7r,) is the vector of residuals measured for the fitted track /. Vi is the covariance 
matrix of the residual measurements of track /. The generic solution for alignment corrections (5a) 
is: 

-1 

5a- 




\i £ tracks 

(3.3) 

where A^' is the covariance matrix for 5a. The size and contents of the matrix A depend on the 
details of the alignment method which are explained in the following sections. 

3.1 The Robust approach 

The Robust alignment approach [6] is an iterative method to align a silicon detector with overlap- 
ping modules. In each iteration alignment corrections are calculated from measurements of mean 
residuals, res, and mean overlap residuals, ovres, in the x and y coordinates. Overlap residuals are 
defined as the difference between two residuals from two overlapping modules, y SCT residuals 
are constructed using both hits from each side in a module. The algorithm only corrects for shifts 
in the plane of the module. The alignment corrections are given by: 



1 



^1 to s-i are defined as: s\ = res;s2 = Y.ovreSx',ST, = Y^ovresy, where 5s j are the measurement 
uncertainties. The range of the sum depends on the geometry of the detector. Given the simple 



CTB geometry, a straightforward implementation of Eqn. 3.4 was used. 
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The alignment corrections were obtained as follows: The fractions of non-overlap and overlap 
hits in the sample were controlled by coefficients A for overlap hits and B for non-overlap hits, to 
adjust the influence of each information on the x and y correction. The corrections were weighted 
with the ratio of the total number of overlap hits A • noh^/y and the number of hits B ■ nh^jy to the 
total number of hits, N^jy. The total residual weight rw^^jy and the total overlap residual weight 
orwx/y obtained this way corresponded to \/5sj^: 

B ■ nhj,/y A ■ nohx/y 
rw,iy = ———, orw,,y = —- (3.5) 

There was one overlap for each two modules in a layer. Thus, this information could be used for 
only one sector which was arbitrarily chosen to be sector 1. The alignment corrections for modules 



in Sector 1 is given by Eqn. and in Sector is given by Eqn. 3.7, as: 



ci-x/y = -orw^iY ■ ovres^/y - rw_,/j, ■ res^/y (3.6) 

ax/y = -r^/y (3-7) 

The CTB alignment was carried out using "unbiased" residuals, i.e. , the hit of the aligned wafer 
on the side of the module was removed from the track fit. About 72,000 events from run 2102355 
were used for the alignment. This run contained about 10 to 50 times more hits than overlap hits. 
Information from residual distributions and overlap residual distributions were weighted so that 
overlap residuals had almost similar influence: setting A to 10 and B to 1 was found optimal. 
Further tests showed that other values affected the speed of convergence rather than the final result. 

There were two major limitations in the application of the Robust algorithm to the CTB data. 
First, significant tilts arose from the hand-mounted modules in the setup. In contrast with the other 
algorithms, the Robust algorithm does not correct for rotations. Therefore, after alignment, the 
residuals still had a global Y dependence, in agreement with the tilts observed around the Pixel y- 
axis (see Sec. The dependence vanished when the modules were rotated accordingly. This is 
the main reason why the residual resolution after the Robust alignment were not as good as the ones 
achieved by other algorithms. The modules with the largest residuals after the Robust Alignment 
correspond to the modules with the largest rotations. Second, discrete Pixel y {rj) residuals resulted 
in less stable mean of the residuals with respect to any small shifts. 

The Robust algorithm converged on a solution without a tight track selection. Although 30 
iterations were performed to align the detector, stable results were achieved after 15 iterations. The 
residuals improved significantly and the track quality stabilized after a small number of iterations. 
After 30 iterations, about 1 jLtm global shifts of module positions in the negative x direction were 
observed (Fig. ^. The Robust algorithm had the advantage of requiring minimal computing re- 
sources. The CPU time used by the algorithm were shown to be negligible compared to that of the 
preceding track reconstruction. 

3.2 The Valencia approach 

The Valencia alignment algorithm [4] is based on the numerical minimization of the x'^ function 



defined in Eqn. 3.2 using "biased" residuals {i.e., the hit of the module being aligned is included in 



the track fit). The covariance matrix is assumed to be diagonal and the diagonal elements are filled 
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Figure 4. Shifts in the x direction for the 30th iteration of the Robust alignment. Statistical uncertainties 
on the measurements are represented by the error bars. The quality of the reconstructed track parameters is 
insensitive to this negligible systematic effect. 



with the measurement uncertainties, a,-, , for residuals, r,, both of which are calculated numerically. 
The only fit parameters are the alignment corrections, ignoring the correlation between track and 
aUgnment parameters. The algorithm is therefore executed iteratively, alternating between track 
and alignment fits. 

The SCT endcap outer module strips follow a fan-out geometry and thus have a variable pitch 
along the vertical direction (Sec. ^. Therefore, instead of using the standard "linear" residual 
(perpendicular distance from the track prediction to the strip), "angular" residuals {8y) were used 
(Fig. H). These represent the difference between the angular separation of the signal channel and 
a "fictitious" strip passing through the extrapolated point. The strip-pitch dependence was thus 
avoided, and uniform angular residuals were obtained. 

The outlier hit rejection was applied by defining an acceptance region determined by a critical 
value of the (outlier rejection). This value was taken as three standard deviations with respect 
to the mean value of the reduced residual distribution (r,/ar, ) calculated before the minimization. 
The fraction of measurements lying out of the acceptance region was ~3%, and reduced to below 
1 % if five standard deviations were used. 

The Valencia algorithm was intended for runs without magnetic field yielding straight tracks. 
After reconstruction, each track was extrapolated to the silicon modules. If the extrapolation lay 
outside the module geometrical acceptance, the track prediction was discarded. The module inter- 
section point of the accepted tracks was transformed into the local frame and residuals were cal- 
culated. For Pixels, only measurements in the {x) direction were considered. The Tj -coordinate 
was ignored due to the non-Gaussian residual distributions (Fig. ||). For the SCT modules, angu- 
lar residuals and measurements from both SCT sides were used^ . Although an analytical residual 
linearisation as a function of the alignment parameters was not computed, the dependence of the 
on the alignment parameters remained linear. Fig. |6| shows the contour regions for two fitted 
variables and three different confidence level intervals (68%, 90% and 95%) for one Pixel module. 

^Except for module [layer 2, phi 1] with a single working plane. 
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Figure 5. Residual definitions in a SCT end-cap module used by the different alignment algorithms: the 
angular and linear residuals and distance of closest approach (DOCA). 
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Figure 6. Contour ellipses of constant probability content for three different confidence levels and two 
alignment parameters of a Pixel module; translation along the most sensitive coordinate (horizontal axis) 
and tilt angle (vertical axis), as obtained by the Valencia alignment algorithm. 



The alignment was performed in three consecutive steps, with variable number of iterations 
in each step: (st-l) internal alignment of the Pixel modules, (st-2) broad alignment of the SCT 
modules with respect to the Pixel system, and {st-3) fine alignment of all silicon modules. In st-l 
(~ 6 iterations), tracks lying in the overlap region between Pixel modules in the same layer were 
selected to enhance the number of overlap hits and to produce a pixel aligimient. In st-2 (~ 2 
iterations), tracks reconstructed only with the pixel hits were extrapolated to the SCT planes. In 
this manner, it was possible to compute SCT residuals (unbiased only in this case) which served as 
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input for an initial alignment of the SCT modules with respect to the Pixel modules. The required 
correction of the SCT modules was several hundreds of microns. In st-3 (~ 8 iterations), all silicon 
modules were included in the track fit and all were aligned simultaneously. In this last stage the 
alignment corrections per module were of few micrometers. 

During alignment, the first Pixel module [layer 0, phi 0] was kept as an anchor; fixed to its 
nominal position to fix the global degrees of freedom. DoFs to which the sensitivity was very small 
were excluded from the set of fitted alignment parameters. Module positions along the beam axis 
were not considered. For Pixels, only the displacements along the sensitive coordinate were fitted. 
The tilt angle (7) was excluded in {st-1 and st-2), but fitted in step {st-3). For the SCT modules, the 
parameters for displacements along and across the sensitive coordinate together with the in-plane 
rotation were fitted in all steps. The inclusion of one additional angle (j8) during the last iterations 
was found to marginally help to improve the results for both sub-detectors. 

3.3 The Local approach 



The Local x approach [7, 8] derives from Eqn. 3.2. The ;^ -function uses unbiased residuals. 



which are defined as the 3D distance of closest approach (Fig. g). The algorithm uses a diago- 
nal covariance matrix, V that is simlar to that of the Valencia approach. The residual errors are 
calculated using hit errors and the extrapolated tracking errors. 

The Local algorithm produces alignment constants for each module separately, neglecting 
correlations between the modules during an iteration. Thus, the solution reduces to inverting as 
many N xN matrices as there are modules, where N corresponds to the DoFs of each module (up 
to 6). Track parameters with a better fit quality gradually bring correlations into play after every 
iteration. 

The fact that CTB was found to be a degenerate setup for track-based alignment required 
inclusion of external constraints to resolve some of the degeneracies. These were a momentum 
constraint to the reconstructed tracks and an additional stabilization term to the diagonal elements 



of the matrix A in Eqn. |3.3| . The stabilization term acts like an additional measurement with a 

is 



zero residual, full sensitivity in the corresponding degree of freedom (the derivative in Eqn. 3.3 
equal to one) and an uncertainty agtabiiity The uncertainty astabiiity corresponds to the inverse of the 
square root of the added term. These additional stability terms constrain the movement to be within 
^stability The valucs for astabiiity are 10, 10, 100 /im for the Pixel x,y,z coordinates and 100 /im for 
the SCT x,y,z coordinates. For the module rotations the value for astabiiity is set to one mrad. 

The momentum of the incident particles from SPS is known more precisely than the intrinsic 
momentum resolution of the CTB ID setup. Consequently, this information can be used to constrain 
the track curvature. Tracks with different beam energies were used as input, using 10,000 events 
from each pion run listed in Table [T[ The alignement procedure was parallelized where multiple 
jobs with different momentum constraint settings were executed simultaneously. When jobs were 
finished the alignment information was collected and merged. Subsequently a new iteration with a 
new set of parallel jobs was started. 

The usage of overlap hits, a hit lying in the overlap region of two modules on the same layer, 
has a profound impact on alignment by constraining relative positions of both sectors, thus avoiding 
divergences due to lack of external constraints. Residual calculation is also more precise for overlap 
hits than non-overlap hits. Non-overlap hits were rejected for alignment once a defined limit was 
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Figure 7. Flow of alignment parameters of the 6 Pixel modules through the iterations of the Local 
alignment algorithm. 



reached. In this way the number of overlap hits was enriched with respect to the number of non- 
overlap hits. The maximum number of non-overlap hits was set to 400. 

For the alignment the iteration chain was performed 60 times. The flow of the 6 alignment 
parameters of each Pixel module through the iterations is shown in Fig. ^. After 10 iterations, nearly 
all degrees of freedom of all modules converged on stable values. Slower convergence of some 
parameters was due to the imposed stability term. The procedure was stopped after 60 iterations, 
when no significant improvement of track parameters was observed and alignment corrections for 
the sensitive coordinates were at the submicrometer level. 

3.4 The Global approach 

The Global x} algorithm [9, 10] is based on the minimization of the x^ defined as Eqn. 3.2 with 



respect to the alignment parameters. The residuals are defined within the module plane and are 
biased (i.e., the hit of the module being aligned is included in the track fit). They depend on the 
track parameters {%) as well as on the subset of alignment parameters related to the intersected 
module (a): 

dx dr dr dn 
da oa on da 

The method has the advantage of properly treating all correlations between residuals arising 
from common track parameters and Multiple Coulomb Scattering (MCS). Solution requires 
inverting a symmetric matrix of size N x N, where N is the number of DoFs of the problem. For 
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Figure 8. Left: original eigenspectrum of the unconstrained solution (84 DoFs - open squares) overlaid 
with that of the actual solution (42 DoFs - solid triangles) of Global algorithm. Six solid circles on the 
leftmost are the singular modes of the asymptotic freedom of the system. Their actual value is zero up to the 
numerical precision and were fixed to 10^^ for clarity only. Right: actual solution eigenspectrum (42 DoFs 
- soUd triangles) compared to the analogous one without the anchor modules requirement (open squares). 
The five solid circles are the new near-singularities. 



large systems (for instance, the entire ATLAS ID), the solution with accurate numerical precision 
and in a reasonable CPU time could be a challenge [24]. In the CTB case, however, the system 
consisted of just 14 silicon modules. Therefore it was free from such numerical limitations. In- 
trinsic alignment of an unconstrained system always leads to a singular matrix and consequently 
an ill-defined solution. This is best solved by diagonalization of the matrix. The singular modes 
can subsequently be ignored in the solution. The procedure can be further extended to remove all 
"weak modes" which either represent unphysical deformations or have an associated error exceed- 
ing expected misalignments. 

In order to solve the CTB alignment, the following approach was adopted: two anchor mod- 
ules were chosen (the first Pixel and the last SCT) which removed the exact singularities from 
the solution. All considered tracks were nearly parallel to one another and orthogonal to the SCT 
module planes. Also the y tilt angles of the Pixel modules were considered to be very accurately 
known from the survey. Consequently the following DoFs were removed from the fit: out of the 
plane translation and the rotations with respect to x and j-axes. This choice resulted in 3 DoFs 
per module (36 in total). However, results indicated a substantial residual misalignment related to 
the uncorrected y rotation of the Pixel modules. The largest misalignment was found for the upper 
module in layer 2 with a value of 25.2 ± 0.5 mrad. The y rotations of the Pixel modules were even- 
tually included in the alignment fit which efficiently eliminated the corresponding misalignments. 

The final alignment was concerned with 42 DoFs. Fig. |8| shows the comparison of eigenspectra 
(obtained by DSPEV routme. from the LAPACK library [25]) of the unconstrained CTB geometry 
and the one used for the final alignment. Elimination of unphysical parameters efficiently removed 
the lowest part of the eigenspectrum. Fig. ^ gives also the comparison of the final alignment to the 
one without anchor modules. The five weak modes correspond to the approximate"* freedom of two 
global translations and three rotations of the entire setup. 

^ Axes of local reference systems in different modules are not parallel which lifts the perfect translational degeneracy. 
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Figure 9. Residuals of all Pixels (left) and all SCT (right) modules before and after alignment. The width of 
the distribution is consistent with the intrinsic resolution of the modules. 



The method required four iterations for convergence, however a total of seven iterations was 
used on about 50,000 events at each iteration. Translations of some SCT modules in y-direction 
were found to be as large as 1.5 mm; x translations never exceeded 0.4 mm. 

4. Results 

In order to assess the quality of the alignment, one must check the track reconstruction quality and 
physics observables. For this purpose, the alignment corrections were applied to the data detailed 
in Table [T[ 

After aligning the modules, the track finding efficiency increased. For example, for the Robust 
alignment approach, the number of tracks per event was found to stabilize at around 0.95. As 
expected, an average of three hits in the pixels and eight in the SCT (two per module) were found. 
All four alignment approaches produced similar performances, consistent with the simulation. 

A well-aligned setup returns residuals with a mean of zero and a width consistent with the 
intrinsic resolution of the detector and the track fit errors. Fig. ^ shows the biased ;c-residuals 
of all the Pixel and SCT modules for the 100 GeV pion run, for those tracks which had at least 
three pixel and six SCT hits. The width of the distribution after alignment is consistent with the 
intrinsic resolution of Pixels and SCT modules. Fig. ^ shows the mean of Pixel module residuals 
for an example run (20 GeV/c pion run). While simulation residual means are centered around 
zero for all modules, the aligned detector data show fluctuations. From the size of the fluctuations, 
we conclude that the Pixel residuals of all alignment methods agree within 5 jim over the whole 



momentum range. Fig. IC also shows a good agreement between the minimization methods and 
the simulation on the residual resolutions. The Robust method results in a worse residual resolution 
since this method only corrects for alignment shifts in the module plane. Fig. |l^ also reveals a 
dependence of the a of the pixel residuals on the module number. This indicates contributions to 
the resolution from the geometry of the setup in addition to the intrinsic detector resolution. The 
residuals also vary because the track error varies along the track due to MCS, for example. 



Similarly, rotations are free only approximately. 
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Figure 10. Pixel biased -residuals mean (left) and G (right) distributions for the aligned data and the 
simulation, for the 20 GeV pion run with B=1.4 T. 
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Figure 11. SCT biased residuals mean (left) and (7 (right) distributions for the ahgned data and the simula- 
tion, for a 20 GeV pion run with B=1.4 T. A single SCT module has two entries in each plot (e.g. 000 and 
001 correspond, respectively, to the module front and back plane at layer 0, sector 0). The dip at index 210 
in the right plot corresponds to the plane which was not working. 



The residual mean distributions for each SCT plane are shown in Fig. 11. Systematic corre- 
lations in the signs of the means is observed among the alignment algorithm results. Fig. 1 1 also 
shows that the residual resolution of the SCT modules in aligned data (except those reconstructed 
using Robust method alignment corrections) are around 20 /im, which is in good agreement with 
the simulation. 

All track parameters at the perigee (do, zo, and Oq and the momentum) were examined when 
tracks were reconstructed with the alignment corrections from the four algorithms. The values of 
spatial track parameters were not exactly similar for tracks reconstructed with different constants, 
however, they followed consistent trends for the runs studied. The difference can be attributed 
to the insufficiently constrained global degrees of freedom. The residuals and curvature, hence 
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the track fit and the pj, are invariant under rigid body translations and rotations of the whole 
system. 

The reconstructed 00 and values depend on the beam properties as well as the module 
locations provided by the algorithms. Therefore, the measured (^o ^^d % in data were used to tune 
the beam spread in the simulation and to evaluate the interaction length, Xq, upstream of the CTB 
setup. The track 0o and % resolutions improved with increasing momentum, as expected with a 
reduced MCS for more energetic particles. 

The momentum reconstruction provides a very powerful test of the alignment performance. 
Fig. |l^ shows the recovery of the momentum resolution of the 100 GeV pion run after alignment, 
from a highly degraded initial measurement. The momentum measurement does not depend on 
global transformations. Therefore the momenta of the tracks reconstructed with different alignment 
constants ought to agree. Fig. 13 is used to compare the electron and pion momenta resolution 
as a function of the reconstructed momentum obtained from the four alignment methods to the 
simulation. The momenta reconstructed using all algorithms, in particular minimization ones, 
are consistent with the simulation. The Robust method returns shghtly worse results since the 
alignment does not take rotations of the modules into account. 

The reconstructed electron momentum is significantly less than the nominal (set by the beam- 
line), for both data and simulation. The presence of several layers of upstream material can account 
for this effect, because the radiated energy of electrons before they enter the tracking volume was 
not recovered. As pions do not suffer as much from energy loss, their reconstructed momenta agree 
much better with the nominal set by the SPS. 

The convergence of alignment corrections per iteration and the improved residual distributions 
presented are mandatory but not sufficient to ensure the success of the alignment. Unfortunately, 
survey data of the CTB detector setup does not exist, therefore, a comparison with the derived 
alignment sets was not possible. However, a comparison of the position and orientation of each 
detector element derived with the four algorithms served as a means of validation. 

When the alignment constants for the four algorithms were compared, the algorithms were 
observed to provided large corrections (several hundred microns). The chosen ahgnment strategy 
(fixing one or several modules as opposed to leaving the whole system free floating or constraining 
or removing some DoFs from the alignment) has an impact on the solution of the global DoF of 
the system. In order to compare the results of the different algorithms, they need to be globally 
matched. Allowing a global offset for each alignment set was chosen to be the method of find- 
ing a best match of the alignment results. After having subtracted the global offsets between the 
geometries, it was observed that a good agreement between the algorithms for the most sensitive 
coordinates x, y and y was reached. Given the low sensitivity of the alignment procedures to the 
alpha and beta rotations, the agreement between the algorithms for these coordinates was only 
marginally improved. 



5. Summary and Conclusion 

Four independent algorithms were used to successfully align the setup formed by the silicon mod- 
ules of the ATLAS Inner Detector tracker, using data collected during the 2004 Combined Test 
Beam. The reconstructed track parameters and hit residual distributions were studied. The per- 



-16- 



CO 
-t— ' 

1 400 

2 350 
I 300 

250 
200 
150 
100 
50 
0, 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



T 



T 



Before alignment 
After alignment 



<^> = 0.0101 (GeV"V 
a(^) = 0.0019 (GeV y 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 

9(GeV-^) 



Figure 12. Momentum resolution of Pixel and SCT detectors for a 100 GeV pion run with non-zero B-field 
with and without alignment corrections. 
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Figure 13. Momentum resolution as a function of the reconstructed momentum. 



formance of the alignment algorithms was assessed by comparing with a simulation, in which all 
modules were at their nominal positions. The simulation can be taken as a benchmark where all 
errors were regarded as only being due to the intrinsic resolution of the modules. 

All alignment approaches yielded results for the reconstructed momentum of electrons and 
pions that agreed with the simulation. Slightly worse momentum resolution was observed using 
the Robust algorithm. This was understood and explained by the fact that the algorithm was lim- 
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ited to re-alignment of the two in-plane translations only. The unresolved residual misalignments 
(e.g. in-plane rotations) unavoidably led to reduced track fit quality and consequently increased un- 
certainties on the reconstructed curvature. For the remaining track perigee parameters, consistent 
results were obtained with each method. 

All four methods agree well on the residuals for all modules and planes, and with the simu- 
lation. The resolution of individual pixel modules is around 10 jxra and the SCT around 20 jxra. 
Observed differences for the residual mean values remain below 5 /xm. We conclude that the sil- 
icon modules of the ATLAS ID were aligned at the CTB with a precision of 5 /im in their most 
precise coordinate. 

The data collected at the ATLAS Combined Test Beam in 2004 served as an invaluable test bed 
for the Inner Detector alignment algorithms. For the first time ever, the readiness of the aUgnment 

algorithms was assessed with experimental data. All algorithms performed satisfactorily given the 
limitations inherent to the CTB geometry and the beamline arrangement. The narrow tower of 
modules and almost parallel particle beams gave rise to undetermined degrees of freedom. These 
were successfully dealt with by the four algorithms, each in its own way, providing consistent and 
high quality measurements of the test beam track parameters. 
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